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ABSTRACT 

We have investigated the report of Tamura and Takeuti that the Ha Hne of the F-type 
post-AGB star HD 56126 is variable on time scales of minutes. To this end, HD 56126 
was observed on two occasions with the William Hcrschel Telescope. Seventeen, re- 
spectively thirty spectra were taken within time span of 1.5 hours in order to detect 
any short term variations. We find that the Ha line profile changed strongly over the 
two month interval, but no evidence is found for short term variability. The variability 
Tamura and Takeuti claim to find is probably due to the low signal-to-noise in their 
spectra. 

Keyvifords: stars: AGE and Post-AGB - stars: circumstellar matter - stars: evolution 
of - stars: mass-loss - stars: individual: HD 56126 



1 INTRODUCTION 



HD 56126 (F5I, 7''13™25.3" 4-10°05'09") belongs to the 
class of the so-called high latitude supergiants that are 
thought to be in the post- Asymptotic Giant Branch (post- 
AGB) phase of their evolution. These are stars with super- 
giant spectral types, but are located at unexpectedly high 
Galactic latitudes. After it was found that most of these 
stars show far-infrared excess emission due to thermally re- 
radiating dust (Parthasarathy and Pottasch 1986; Trams et 
al. 1991; Oudmaijer et al. 1992; Kwok 1993) the post-AGB 
nature became a widely accepted interpretation for these 
stars. The infrared excess is explained as being a remnant 
of the high mass-loss the stars underwent during the AGB 
phase. The evolved status makes the objects less luminous 
and older compared to massive supergiants, placing them 
closer to the plane and allowing them more time to get there, 
solving the high latitude problem. 

During the AGB phase, carbon rich material from the 
helium burning shell was dredged up to the photosphere of 
HD 56126. This is suggested by the detection of the 3.3 fim, 
infrared "PAH" feature by Kwok et al. (1990), and the de- 
tection of C2 and CN absorption lines in the optical spec- 
trum by Bakker et al. (1995). A study by Parthasarathy et 
al. (1992) showed that HD 56126 is metal deficient by a 
factor of ten with respect to solar metalicities, although C, 
N, O and S appeared solar. This photospheric abundance 
pattern resembles that of interstellar gas, where the metals 
have been condensed onto grains (Lambert et al. 1988; Bond 
1991; Van Winckel et al. 1992). 

One of the main problems in the study of late stages of 



stellar evolution of low to intermediate-mass stars is whether 
there is mass-loss during the post-AGB phase, and if there 
is, the magnitude of the mass-loss rate. After the core mass, 
the post-AGB mass-loss rate is the most important parame- 
ter governing the time scale of the evolution of these objects 
from the AGB to the Planetary Nebula (PN) phase. For ex- 
ample the transition of a star with a core mass 0.546 M© , 
having a Reimers mass- loss of 10~* M0 yr~^ takes 100,000 
years (Schonberner 1981 and 1983). A period in which the 
AGB wind has dispersed into the interstellar medium long 
before the central star has reached a temperature high 
enough to ionize its circumstellar shell and become observ- 
able as a PN. However, as pointed out by Trams et al. (1989), 
increasing the mass-loss rate to 5 x 10~^ yr~^ accel- 

erates this transition to only 5000 years, which would make 
a PN readily observable. It is difficult however to determine 
post-AGB mass-loss rates. All "classic" tracers of mass-loss 
rates such as circumstellar CO millimeter emission, far in- 
frared excesses etc. trace the former AGB wind rather than 
the present-day mass-loss in post-AGB stars. 

Most of the post-AGB stars exhibit mild Ha emission. 
Some objects show strong P-Cygni type emission, while most 
of the stars have shell type profiles. In the recent literature 
the Ha emission has been interpreted as the star undergoing 
a post-AGB mass-loss (Trams et al. 1989; Parthasarathy et 
al. 1992; Slijkhuis 1992). In fact, the Ha emission is, next 
to the near-infrared excess observed in some of the objects 
(Trams et al. 1989), the only diagnostic suggested for post- 
AGB mass-loss. Although a P-Cygni type emission can be 
understood with a simple expanding wind model, the dou- 
ble peaked shell type profiles require somewhat more com- 
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plicated geometries as for example rotating disks (Waters et 
al. 1993). The mass- loss interpretation gets more into dif- 
ficulties when considering the fact that many of the stars 
show variable Ha profiles. Some of the stars have shell type 
Ha lines where the ratio between the blue and the red 
peaks varies on time scales of months (Arellano Ferro 1985; 
Waelkens et al. 1991, Waters et al. 1993). Only in the case 
of HR 4049 the Ha variability could be correlated with ra- 
dial velocity variations and the orbital phase of the star in 
a binary system (Waelkens et al. 1991). Although the num- 
ber of binary stars in the sample of post-AGB stars under 
consideration is relatively high (Waelkens and Waters 1993), 
not all stars that exhibit Ha variations are yet identified as 
members of a binary system. 

Another interpretation that can give rise to variable 
Ha emission of post-AGB stars stems from stellar pulsa- 
tion theory. The F-G type post-AGB stars are located in 
a part of the HR diagram where the photosphere is sub- 
ject to instabilities (Sasselov 1993; Zalewski 1992). A conse- 
quence of such instabilities has been pointed out by Lebre 
and Gillet (1991a-|-b, 1992) in their work on RV-Tauri and 
W-Virginis stars, objects that resemble post-AGB stars very 
closely. These stars show variable Ha emission during their 
pulsation cycle. Lebre and Gillet showed that inward and 
outward moving shocks within the photosphere can produce 
line profiles remarkably similar to those observed in the post- 
AGB stars discussed above. 

In this work we aim at confirming the rapid changes 
of the Ha profile of HD 56126 described by Tamura and 
Takeuti (1993, hereafter TT93). These authors reported 
Ha line profile variations on time scales of the order of 
minutes. If the variability on minute time scales is real, it 
will prove very hard to attribute the observed Ha emission 
to mass-loss. In principle it may be possible to have short 
"puffs" of mass ejected from the photosphere and falling 
back again, giving rise to blue and red emission peaks re- 
spectively. Such "puffs" are then indicative of a turbulent 
photosphere, so the mass-loss and stellar pulsation should 
be closely intertwined. It would then be more reasonable to 
explain the rapid Ha variability in terms of pulsation solely 
than a combination of mass-loss and pulsation together. Re- 
alizing the important implications caused by the short term 
variability reported by TT93, we decided to obtain a set of 
high quality spectra with short exposure times of HD 56126. 
For this purpose the Utrecht Echelle Spectrograph mounted 
on the 4.2m WHT telescope on La Palma was employed. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

The first set of observations was carried out on December 
21"* 1993 in service time with the Utrecht Echelle Spectro- 
graph (Unger 1994) mounted on the Nasmyth platform of 
the 4.2 m William Herschel Telescope, La Palma, Spain. 
The weather was fair, but the seeing was larger than 2", 
which, compared with the slit width of 1", resulted in a 
somewhat lower signal-to-noise of the spectra than could 
have been possible. The detector was a 1124 x 1124 TEK 
CCD. Wavelength calibration was performed by observing a 
Thorium- Argon lamp. The central wavelength of the setting 



Table 1. Log of the observations December 21 1993 



Time tcxp Airmass SNR 

U.T. [scc.l 



23:24:39 


60 


1 


.465 


50 


23:27:58 


60 


1 


.446 


52 


23:30:56 


90 


1 


.430 


62 


23:37:41 


90 


1 


.394 


63 


23:41:05 


90 


1 


.378 


63 


23:44:40 


180 


1 


.361 


87 


23:49:35 


180 


1 


.339 


89 


23:54:29 


180 


1 


.318 


83 


23:59:23 


180 


1 


.299 


82 


00:05:25 


240 


1 


.276 


98 


00:11:19 


240 


1 


.256 


90 


00:17:10 


240 


1 


.237 


99 


00:23:09 


240 


1 


.219 


95 


00:29:12 


240 


1 


.202 


103 


00:35:08 


240 


1 


.186 


100 


00:41:11 


240 


1 


.171 


90 


00:47:11 


240 


1 


.158 


94 



SNR has been calculated in the line free region between 6533-6540 A 

was 5587 A, resulting in a wavelength coverage from 4700- 
7200 A. In order to make a compromise between signal-to- 
noise and time resolution, it was decided to make exposures 
with increasing integration times. 

A second set of data was obtained during a run on the 
WHT in the night of February 26-27 1994. The settings 
were the same as for the December run, except for the cen- 
tral wavelength (7127 A), resulting in a wavelength coverage 
from 5500 A, - 1 fim. This time the strategy was aimed at 
searching for time variations within minutes. Thirty expo- 
sures of 60 seconds and 30 seconds were taken consecutively. 
The sky was cloudy, resulting in many telluric features in the 
spectra. The spectral resolution as measured from telluric 
absorption lines (full width at half maximum) is 8.5 km s~^ . 
The logs of the observations and the resulting signal-to-noise 
ratios (measured in the Ha order) are provided in Tables ^ 
and ^ 

2.2 Data reduction 

The data were reduced using the standard Echelle reduc- 
tion software in the IRAF package. The data were bias 
subtracted and divided by a normalized flatfield frame. Af- 
ter the fiatfielding, the individual orders were optimally ex- 
tracted from the image frame, where, using simple photon 
statistics, pixels are weighted. 

The spectra contained 44 orders in the December run, 
and 48 in the February run, that were continuum normal- 
ized for further analysis. The wavelength calibration was 
done by fitting a S*^"* order Legendre polynomial in the dis- 
persion direction and a 4"* order polynomial in the cross- 
order direction. The dispersion of the fit through the more 
than 500 identified Th-Ar lines was 3.8 mA and 7 mA for 
the December and February spectra respectively. During the 
process of extraction the signal-to-noise ratio of the pixels 
was calculated with simple photon statistics. Checks with 
the actual observed signal-to-noise of the data proved to be 
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Table 2. Log of the observations February 27 1994 



Time 


tcxp 


Airmass 


SNR 


U.T. 


[sec] 






23:30:35 


60 


1.133 


58 


23:34:37 


60 


1.141 


55 


23:37:39 


60 


1.147 


53 


23:40:39 


60 


1.153 


37 


23:43:52 


60 


1.160 


47 


00:22:18 


60 


1.26 


72 


00:25:39 


60 


1.280 


60 


00:28:51 


60 


1.292 


68 


00:31:54 


60 


1.303 


55 


00:35:23 


60 


1.317 


60 


00:38:37 


60 


1.331 


48 


00:41:48 


60 


1.344 


55 


00:44:58 


60 


1.359 


73 


00:47:53 


60 


1.372 


52 


00:56:58 


60 


1.417 


64 


01:00:10 


60 


1.435 


57 


01:03:23 


60 


1.453 


60 


01:06:38 


60 


1.471 


58 


01:10:06 


60 


1.492 


53 


01:13:26 


60 


1.513 


59 


01:16:39 


60 


1.534 


50 


01:19:57 


60 


1.557 


60 


01:23:08 


30 


1.580 


45 


01:25:32 


30 


1.597 


29 


01:27:57 


30 


1.616 


33 


01:30:39 


60 


1.637 


54 


01:33:32 


60 


1.660 


63 


01:36:26 


60 


1.685 


56 


01:39:21 


60 


1.711 


59 


01:42:18 


60 


1.738 


53 



correct within 15%. The SNR spectra calculated in this way 
will be used further in this study. 



3 IS THERE VARIABILITY? 
3.1 Time scales of minutes 

In Fig. |l] all individual spectra around Ha obtained Decem- 
ber 1993 are plotted. In order to decide whether there are 
any variations in the spectra we adopt the simple, yet pow- 
erful statistical formalism presented by FuUerton (1990) and 
Henrichs et al. (1994). With this method the variability can 
be expressed in a temporal variance spectrum {TVS, Eq.|l|): 



{TVS)> 



N 



N 



MX) 



(1) 



where A'^ is the number of spectra, Fav(A) represents the 
constructed average spectrum, Fi{X) the individual spectra, 
and ai(X) = Fi{X)/{SNR) of each individual pixel of the 
spectra. 

Following Henrichs et al. (1994) and FuUerton (1990) 
the temporal sigma spectrum ( TSS = VTVS) is calculated. 
This quantity represents approximately (aoba/cTnv), that is 
the standard deviation of the variations of the individual 
spectra with respect to the average spectrum divided by the 




6B76 



6877 



6878 



6879 



A (A) 



Figure 2. An overplot of the first spectrum (thick line, airmass 
1.47) and the last spectrum (thin line, airmass 1.16) taken around 
6875 A in December is shown in the upper panel. The lower 
panel shows the Temporal Sigma Spectrum. It can be seen that 
the variable telluric absorption lines show up as significant peaks 
in the TSS 



standard deviation of the average spectrum. If no significant 
variations are present in the spectra, the value will be close 
to one, significant deviations are directly represented in units 
of the noise level, that is to say a peak "Temporal Sigma 
Spectrum" of three corresponds to a variability at a 3a level. 

The individual spectra and the resulting TSS are plot- 
ted around Ha in Fig. 0. It is clear from the temporal sigma 
spectrum that no significant variability is present. The spikes 
in the TSS spectrum correspond to cosmic ray hits that re- 
sulted in higher count levels and consequently higher signal- 
to-noise ratios. These cosmic ray events are still a bit present 
in the spectra, and provide an independent check whether 
the method is indeed capable of detecting variability. An- 
other way to demonstrate the adequacy of the method can 
be found in Fig. |2|, where we have plotted the first and last 
taken spectrum around 6875 A during the December run. 
The effect of the decreasing airmass on the strength of tel- 
luric absorption lines is seen. 

The same procedure was carried out for the February 
observations. The average spectrum and the resulting TSS 
are shown in Fig. ^. No variability can be found either in 
the February run. The shallow shape of the TSS arises from 
the fact that the normalization of the spectrum was difficult, 
leading to a small variation of the continuum in the wings 
of the Ha line. 

Similar calculations were performed for spectra with the 
same integration times, and for the rest of the UES spec- 
trum. Except for strong telluric absorption lines (see Fig. ^ 
no variations are present. 
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6540 6550 6560 6570 8580 



A (A) 



Figure 1. The 17 spectra that were obtained around Ha during the December run. The spectra are continuum normaUzed, each spectrum 
is plotted with an offset of 0.3 for comparison. The Temporal Sigma Spectrum (see text for details) is also plotted. Note that the small 
spikes in the individual spectra are immediately visible in the TSS 




6540 6550 6560 6570 6580 



A (A) 

Figure 3. The summed spectrum and the resultant Temporal 
Sigma Spectrum of the observations obtained in February 1994. 
The slight curvature of the TSS arises from the fact that the 
normalization of the spectrum was difficult, leading to small vari- 
ations of the continuum in the wings of the Ha line 



3.2 Time scales of months 

Average spectra were computed by adding all spectra be- 
fore the continuum correction, resulting in spectra with a 
high SNR{> 300 and < 175 for the December and Febru- 
ary observations respectively) . The radial velocity is not sig- 



nificantly variable, the heliocentric velocity we derive from 
strong metallic lines are 87±2 km s"^ , and 84±2 km s"^ 
for the December and February spectra respectively. 

In Fig. ^the Ha lines at both occasions are shown. The 
profiles are different; whereas in the December spectrum the 
Ha profile shows a central absorption at 103 km s"'^ , it 
has shifted to 94 km s~^ in February, the peaks around the 
Ha absorption have increased in strength. The true shape of 
the emission is hard to determine, but considering the red- 
shifted center of the Ha absorption, there appears to be a 
rather broad emission component blue-shifted by more than 
20 km s"^ . 

Not only the Ha lines differ, a preliminary comparison 
of the individual spectra shows that all spectral lines are 
broader by more than 40% at their full width half maximum, 
while the depression relative to the continuum decreased in 
the February spectra relative to the December spectrum. 
Interestingly, except for the strongest metallic lines (depres- 
sion < 50%) where a blue-shifted absorption wing is visible, 
all lines appear symmetric. The interstellar component in 
the profiles of the Nal Dl & D2 lines do not change in ve- 
locity nor in shape between the two spectra, indicating that 
changes in other lines are not due to observational errors. 
For the Fel and Fell lines we note that the equivalent width 
of the Fel lines decreased, while the equivalent width of the 
Fell lines has increased. The fact that the ionization de- 
gree of Fe has changed indicates an increase of the effective 
temperature, the broader lines are probably the result of a 
higher surface gravity. The deeper central absorption in the 
Ha line is also due to an increase of the effective tempera- 
ture which gives a higher population of the n = 2 level for 
hydrogen according to the Boltzmann law. 
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Figure 4. Ha profiles of HD 56126 at the two epochs. The dif- 
ference spectrum is shown in the lower panel 



loss, where the stronger emission in the February spectrum 
is simply due to the increased temperature of the star, it 
could also be due to a shock wave ploughing through the 
photosphere. It may be worthwhile to monitor post-AGB 
stars on time scales of days and weeks to appreciate the 
true temporal changes in their spectra. 

Finally we consider whether the absence of any short 
term variability can be expected based on simple arguments; 
in principle the fastest time scale in which a wave can prop- 
agate in a stellar atmosphere is the scale height h divided 
by the speed of sound Vs : 



h ■ 



kTW 



kT 



[km s 



(2) 



(3) 



With M the mass of the star and R its radius, fi is the mean 
molecular weight which is taken to be 1. Adopting typical 
parameters for a F-type post-AGB star, M — 0.6 Mq , and 
R — 6O7?0, we obtain a shortest possible time scale of 1.8 
days. 

One can also estimate the pulsation period of a typical 
post-AGB star: 



4 DISCUSSION 



rp 



In Sect. 3.1 it was shown that no short term variability is 
found in the Ha line of HD 56126. 

Before we discuss any explanation for this observa- 
tional result, let us first review the data presented by TT93. 
They presented five spectra of HD 56126 taken with the 74 
inch Okayama telescope. The spectra were obtained within 
47 minutes, with 6 minutes the shortest time interval be- 
tween two observations. The wavelength resolution is about 
0.13 A (6 km s~^ at Ha ), comparable to the spectra pre- 
sented here. From their plots we find for their best spectrum 
a peak-to-peak noise of approximately 17% of the continuum 
level and of about 30% for the other spectra. We estimate 
their SNR at about 10-20. Of the five spectra three show 
an absorption profile, only two show the variation noticed 
by TT93. One of these shows a red emission peak while 
another spectrum shows a blue peak, which they claim to 
be "drastic changes of the Ha profile in the order of ten 
minutes". In fact, these are the only spectra that show no- 
ticeable variations. The signal-to-noise of their spectra is at 
most 20, so that a la variation has to be at least a 5% 
change for an unresolved line. The (variable) emission com- 
ponents hardly exceed this value. Furthermore, the strong 
Sri line at 6550.2 A (W = 7 mA, depression to 0.8 x con- 
tinuum) is present in only two of their spectra. This further 
suggests that the observed changes may be artifacts of the 
noise rather than true changes in the spectrum. 

In Sect. it was argued that, due to the broadening of 
most stellar lines and the change in ionization degree of Fe, 
the star has become hotter over a time interval of 65 days. 
This indicates that the star is variable, consistent with the 
findings of Bogaert (1994), who found photometric variabil- 
ity of HD 56126. 

It is still hard to identify the mechanism responsible 
for the emission observed in Ha . It can be due to mass- 




(4) 



Where Q is the pulsation constant in days. Adopting the 
same stellar parameters as above, and noting that Q ranges 
from 0.04 for classical Cepheids to 0.16 for W-Virginis stars 
(Allen 1973) we obtain an expected period of 30 - 96 days. 
This simple exercise tells us that variations on time scales 
shorter than days are not really to be expected in post-AGB 
stars. It is consistent with the change observed over a period 
of two months. 



5 CONCLUSIONS 

We observed the post-AGB star HD 56126 (F5I) at short 
time intervals in order to look for very short term variations 
of the Ha line. Despite the high quality of the data we were 
not able to find significant variations. The report by Tamura 
and Takeuti (1993) that the line is variable within minutes 
could not be confirmed. We believe that the variability they 
claim is due to the low signal-to-noise ratio of their spectra 
and not due to real variability. Strong variations are found 
over the two month interval, where the central absorption of 
the Ha line has shifted and both peaks around the central 
absorption have increased considerably. It is shown that the 
star has increased its surface temperature over two months, 
indicating that the underlying star rather than the circum- 
stellar environment changes, adding more validity to pulsa- 
tional models that account for observed changes in spectra 
of post-AGB stars. 
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